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Abstract 

In this paper we continue the investigation of the effect of quan¬ 
tum number conservations of pions, kaons, and protons, with very high 
transverse momenta (up to 25 GeV/c), during parton fragmentation and 
hadronization in p-p and Pb-Pb collisions at LHC energies ]T][2]. The 
strength of the conservation effects are studied by identified two-particle 
correlations in Monte Carlo generated events in the mid-rapidity region 
(|?j| < 1). The simulated p-p events were generated with PYTHIA 8 [5], 
using its main default settings, at y/s = 200 GeV, y/s = 2.76 TeV, 
yfs = 7 TeV, and y/s = 14 TeV. In parallel to this, HIJING 1.36 [6 
was used to generate Pb-Pb events at ^/snn = 2.76 TeV with centrali¬ 
ties 0 — 10%, 30 — 40% and 80 — 90%. The extracted identified associated 
hadron spectra for charged pion, kaon, and proton show identified trigger- 
hadron dependent splitting between oppositely charged associated particle 
species in any nucleus-nucleus collisions. The Pb-Pb data exhibits a pe¬ 
culiar splitting pattern as a function of the transverse momentum of the 
associated particle pr, assoc both on the near and away side that is dif¬ 
ferent from the patterns observed in p-p collisions. The splitting shows 
smooth evolution with collision energy and event multiplicity in p-p col¬ 
lisions while in Pb-Pb collisions different trend were observed for kaons 
and protons. 


1 Introduction 

Two-particle angular correlation is an efficient and powerful tool to study the 
properties of the underlying physics processes of particle production in different 
collision systems, like proton-proton or heavy-ion collisions. Using identified 
particles we can check those kind of conservation laws which are related to the 
different quark and flavour content of the processes 0. This can help us to 
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understand the specific hadronization processes in more details. Furthermore, 
focusing on the correlations of identified particles one can study recombination 
effects 0 S] besides the parton fragmentation in the hadronization mechanisms. 
We used in our analysis the most known Monte Carlo generators (PYTHIA [5] 
and HIJING 0) of the high energy physics community. This can be practically 
helpful to verify the validity and the prediction power of the studied physical 
quantities of this study to experimental data. 

In this paper we continue our earlier work [2] on the investigation of the 
particle production and the quantum number conservation of identified hadron- 
triggered angular correlation of identified hadrons at high transverse momenta, 
up to pt = 25 GeV/c. This time the strength of the conservation effects is 
studied as a function of collision energy and event multiplicity in p-p collisions 
and as a function of centrality in Pb-Pb collisions. 

This proceeding is organized as follows: Sec. [2] briefly presents the Monte 
Carlo data samples and analysis settings which were used in this study. In Sec. [3] 
we show the obtained results for relative yields of unlike-sign associated particle 
pairs in specific physical configurations in p-p and Pb-Pb collisions. Finally, 
Sec. 2] summarizes the results of this proceeding and gives outlook for possible 
further analysis of the subject. 

2 Monte Carlo simulations of identified 
two-particle azimuthal correlations in 
p-p and Pb-Pb collisions 

2.1 Monte Carlo data sets and analysis cuts 

As in our earlier study, here we also consider generated identified particles: 7r ± , 
K ± , p, p, and the charged hadrons (h ± ) as a reference [2]. 

Event generators PYTHIA 8 0 and HIJING 1.36 [6j were used for the pro¬ 
duction of p-p and Pb-Pb event samples. PYTHIA 8 was used with its default 
settings (tune 4C), while HIJING event generation included the simulated ef¬ 
fects of quenching and shadowing. The detailed data sets for the generation of 
p-p and Pb-Pb samples can be found in Table [Q All charged final-state parti¬ 
cles were kept for analysis from the Monte Carlo events, nevertheless the trigger 
particles and associated particles used for the two-particle correlations measure¬ 
ments were limited in pseudorapidity in order to take into account the sensitive 
geometry of today’s experimental detectors. Trigger particles with \ptrig\ > 0.5 
and the associated particles with |r/ assoc | > 1 were rejected. In order to exclude 
the low-momentum region from the bulk effects, the transverse momenta of the 
trigger particles were considered above pr,trig > 2 GeV/c in all the cases. The 
associated particle momenta were chosen such that pt, assoc < PT,trig to avoid 
double counting. The uncorrelated background was subtracted by applying the 
ZYAM method 0. 

2.2 The definition of the associated per trigger yield 

Measuring quantum number conservation requires constraints both on the trig¬ 
ger and the associated particle directions to tighten the sensitivity of the mea- 
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Collision system 

y/SNN 

Statistics 

P-P 

14 TeV 

500 M 

P-P 

7 TeV 

100 M 

P-P 

2.76 TeV 

100 M 

P-P 

200 GeV 

500 M 

Pb-Pb (0-10%) 

2.76 TeV 

4 M 

Pb-Pb (30-40%) 

2.76 TeV 

2 M 

Pb-Pb (80-90%) 

2.76 TeV 

10 M 


Table 1: The list of Monte Carlo generated data sets used for our analysis. 
Proton-proton events were generated with different collision energies, while Pb- 
Pb events simulated with different centralities as indicated in the parenthesis. 
See the text for further details. 


surement. The associated per trigger yield is defined with the following quantity 
as a function of pseudorapidity difference Ar) = ptrig — Vassoc and azimuthal an¬ 
gle difference Acj) = (j> trig - 4> a ssoc- 

d A" (\ m a \_ 1 d-ZV assoc ( , 1 ^ 

d(Ajj)d(A0) &V,PT,trig, PT,assoc) ~ ^ ' d(A?/)d(A(/)) > W 

where IVtrig and iV assoc is the number of trigger and associated particles per 
event, respectively. This quantity was measured in several pr,trig and pt, assoc 
intervals. The azimuthal angle correlations in this study are always projected 
within the pseudorapidity difference | Ary| < 1. The associated per trigger yield 
on the near side was extracted from the A <f> interval |A</>| < tt/2 and the away 
side from the A0 interval 7r/2 < A <f> < 3tt/2. 

2.3 The PID-associated spectra 

To test the quantum number conservation the identified associated spectra of the 
associated particles (hereafter called PID-associated spectra) have been plotted. 
We have investigated the PID-associated raw spectra using K ± 1 p, p, and 
charged hadron triggers. 

Figure |T] shows the identified particle yields dN^/ 3 ^ c /dpr,assoc for the iden¬ 
tified associated particles up to pt, assoc < 26 GeV/c in p-p collisions at y/s = 
14 TeV. The kaon (green symbols) and proton (red symbols) spectra are scaled 
for better visibility. The PID-associated spectra have been plotted for positively 
charged pions, as trigger particles, selected from the transverse momentum range 
2 GeV/c < p T ,trig < 25 GeV/c. 

On the left side of Figurc[l]one can see the pt, assoc spectra for pions, kaons 
and protons for the ’near side’ case. Right side shows the spectra for the ’away 
side’ case. Full symbols refer to particles with positive charge and negatively 
charged particles are indicated with open symbols. 

The PID-associated spectra of the different associated particle species show 
the expected evolution with the transverse momentum. The yields of the PID- 
associated spectra significantly decrease with the selection of charged pion, kaon 
and proton triggers. 
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Figure 1: (Color online) 7r + -triggered associated particle spectra for near- and 
away side particles scaled by integral of the charged hadron identified associated 
spectra in simulated p-p events at y/s = 14 TeV. Additional scaling was applied 
to the spectra for better visibility. See the text for the details. 


Trigger / Assoc 

Strength of effect 


Near side 

Away side 

7T~*” / (7T+, 7T ) 

~ 2 

~ 1 

K+ / (K+,K~) 

~ 5 

~ 1 

P / (p,P) 

~ 10 

~ 1 


Table 2: Qualitative yield differences between positive and negative associated 
particles on the near- and away sides, listed here on Fig. [I] 


On a qualitative basis one can see on the left side of Fig. |Tj the yield dif¬ 
ferences between the positive and negative associated particles. Looking at the 
right side of the figure, no clear difference can be observed between associated 
particles on the away side. We listed these qualitative observations in Tabic [2] 
It shows essentially that the largest yield difference we observe is for protons 
and the smallest one is present for pions, as open and full symbols are splitting 
towards higher px values. 

We note here, that similar plots can also be obtained using other type of 
trigger hadron species, as we do in the next subsection 12.41 


2.4 Identified particle ratios 

This part was already discussed in Ref. [2] in detail, thus here we remind the 
reader of the main definition and features of the defined quantity. 

The charge (C) and flavour specific or other quantum numbers, like baryon num¬ 
ber (B) and strangeness (S'), of generated particles are expected to hold during 
the hadronization and its traces can be seen on the final PID-associated spectra. 
In order to quantify and magnify the expected quantum number conservation 
effects, we define the following ratio: 


d N, 


PID 


dN. 


hadron 

assoc 


d PT,a 


d PT,a 


( 2 ) 
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This quantity can be understood as follows: in the nominator we consider the 
case when we are interested in the yield of the identified associated particles 
having an identified particle as trigger. While in the denominator we calculate 
the yield of identified associated particles without any PID information on the 
trigger particles, i.e. using charged hadrons as triggers. 



Figure 2: (Color online) PID-triggered associated particle yields relative to the 
charged hadron-triggered associated yields. The upper row shows the PYTHIA 
simulation in p-p and the lower row shows the 0 — 10% central Pb-Pb HIJING 
events. The left column corresponds to the near side (tt/2 < |A(()|) correla¬ 
tions and the right column corresponds to the away side {n/2 < A tfi < 37r/2) 
correlations. See the text for details. 

In Figure[3]the ratios, defined by Eq. ©, for six identified trigger / identified 
associated particle species are shown: 7r + /7r + , tt + /tt~ , K + /K + , K + /K~ , p/p , 
and p/p. We used only the like-sign (+/—), unlike-sign (+/+) particle combi¬ 
nations, since we are interested in the flavour balance between the associated 
particles regarding specific quantum number. 

The upper panel shows the ratio of the yields defined by Eq. (© in p-p collisions 
for the near- and away side, on the left and right respectively. An obvious split¬ 
ting effect can be visible between the unlike-sign associated particle pairs, which 
is the largest for p/p , and the smallest for tt + /tt~ . In contrast, on the away side 
the splitting effect is not really present. This indicates that in PYTHIA gener¬ 
ated particles does not have connection to the choosen trigger particle in this 
context. From this picture we can qulitativley say that the strength of the 
quantum number conservation effects is increasing in the order of charge (C), 
strangeness ( S ) and baryon number {B) respectively. 

Eq. © is plotted for HIJING generated events in the lower panels of Fig. [5J 
as well. A peculiar pattern is present on the near- and away side in the pt, assoc = 
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2 — 8 GeV/c and pr, assoc = 2 — 10 GeV/c ranges for charged kaons and protons, 
where the baryon/meson anomaly was observed at RHIC and at the LHC [S . An 
opposite trend of the splitting at higher transverse momenta can be attributed 
to the minijet production embedded in the generator. 


2.5 Identified particle ratios in different trigger px bins 

So far the associated particle spectra were studied in case of wide range of 
trigger particles’ momenta 2 GeV/c < pr,tri g < 25 GeV/c. Here we briefly re¬ 
view the pr,trig dependence of the spectra. For this we plot the differences 
between the like-sign and unlike-sign particle pairs for the normalized yields 
(dN^c/^PT,assoc)/{AN^ c on /dp T ,assoc) as plotted in Fig. [2] Since positively 
charged triggers were choosen we can introduce the definition below in order to 
plot these ratios higher than one: 


/ dN PID 

I assoc 

y dpT,assoc 


,4 j\rhadron 
assoc 


dp T ,a 


(+/-) 


/ dN PID 
I assoc 

\ dpT,assoc 


d N pid 

’assoc 


dpr,a 


( + /-) 


j Ajhadron ' 
'- Lly assoc 

(+/+) 

dpT,assoc y 

) 

' dN PID N 

assoc 

(+/+) 

dpT,assoc i 

) 


(3) 


where (+/—) and (+/+) denote the unlike-sign and like-sign trigger/associated 
particle pairs, respectively. This double ratio is plotted for four different pr,trig 
trigger bins on the horizontal axis of Fig. 0 in cases of p-p and most central 
(0-10%) Pb-Pb collisions. Deviation of the value from 1 means the relative 
difference between triggered like- and unlike-sign hadron yields. 

In the upper row of Fig. [3] near side ratios are plotted. One can observe 
how the splits - in agreement with Fig. [2] evolve as a function of pt, assoc in the 
near side cone. In general, the splitting effect is present in all choosen pr,trig 
bins, 3 — 7 GeV/c, 7—10 GeV/c, and 10 — 20 GeV/c, with the same order of 
magnitude, for the specific partilce species. We can summerize the observed 
effects as follows: 


• In p-p collisions the effect is present for pions, kaons, and protons in the 
near side cone and the largest for protons with a factor of 10 effect. 

• In Pb-Pb collisions the effect is present for kaons and protons in the near- 
and away side cones and the largest for protons. 

• In Pb-Pb collisions the observed trend of the effect for kaons and protons is 
completely opposite compared to the situatuion in p-p collisions in the in¬ 
termediate associated transverse momentum spectra 3 GeV/c < pr,assoc < 
7 GeV/c interval. 

In p-p collisions at the away side the effect is increasing with px but the 
magnitude saturates to a constant value. A similar, but larger effect is seen in 
central Pb-Pb collisions within the 3 — 10 GeV/c transverse momentum range, 
where the Cronin effect DJi and the baryon/meson anomaly act. 
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Figure 3: (Color online) Yield differences between the oppositely charged trig¬ 
ger/associated particle pairs. The yield enhancement is relative to the positively 
charged trigger particles. Different colors refer to different Monte Carlo data 
sets as indicated on the plot. Upper panel refers to near side, while lower panel 
refers to away side. 


3 PID-triggered unlike sign-to-like sign associ¬ 
ated spectra ratios in p-p and Pb-Pb collisions 

3.1 Energy dependence in p-p collisions 

In Fig. H]we plotted the unlike-to-like sign relative yield ratio defined by Eq. Q 
for pions, kaons, and protons with different collision energies yfs. At the near 
side (upper panel of the figure) one can see a clear dependence of the splitting 
on y/s as a function of pt, assoc for all particle species. As we go up with yfs the 
effect shows a decreasing trend. This can be related to the fact that the higher 
the collision energy the more the number of hard interactions. This has effect 
on the splitting and it is ’washed out’. 

On the away side we have no such kind of yfs dependence, which can be expected 
based on the argument we had in Sec. [2j 
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Figure 4: (Color online) Yield differences between the oppositely charged trig¬ 
ger/associated particle pairs for different collision energies. 


3.2 Event multiplicity dependence in p-p collisions 



1 10 1 10 1 10 


PT, assoc [GeV/c] 

Figure 5: (Color online) Yield differences between the oppositely charged trig¬ 
ger/associated particle pairs for different event multiplicities in = 14 TeV 
p-p collisions simulated by PYTHIA 8 (tune 4C). 




































































In Fig. [3] one can see the event multiplicity dependence of Eq. J3]) for pions, 
kaons, and protons at y/s = 14 TeV in p-p collisions. We observe the same kind 
of trend in the evolution of the splitting with the multiplicity what have been 
seen in the previous case, since the number of generated particles depends on 
the collision energy y/s. This case can be treated as a subset of the previous one, 
becasue here we fix the y/s and this way the generated number of particles are 
bouded. We can conclude that the higher the event multiplicity the smaller the 
size of the splitting due to increasing underlying event contribution at higher 
multiplicities. 

Very similar to the case with y/s dependence we state here that the effects on 
the away side have no obvious dependence with event multiplicity. 

3.3 Centrality dependence in Pb-Pb collisions 



PT, assoc [GeV/c] 

Figure 6: (Color online) Yield differences between the oppositely charged trig¬ 
ger/associated particle pairs for different centralities in Pb-Pb collisions at 
yJsNN =2.76 TeV. 


In Fig. [6] the centrality dependence of Eq. ([3]) is plotted in Pb-Pb collisions 
at y/SNN = 2.76 TeV. We expected here a similar effect such as obtained in the 
previous subsection due to the fact that the multiplicity is related to centrality. 
In peripheral collisions (black symbols, centrality 80 — 90%) we observe a very 
similar behaviour for pions as we have seen in Subsec. 13.11 for the collision 
energy dependence in p-p collisions. For kaons and protons there are two major 
deviations. Firstly, Eq. m on the near side both for the kaons and protons 
starts at the same value as for pions, while in Fig. [3] at the same collision energy 
(-y/s = 2.76 TeV) we see an obvious inverse deviation from unity in the first 
Pt, assoc bin-which has the highest value for protons. This deviation turns to 
be greater than 1, approaching to highest pt values. Secondly, we do see an 
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inverse evolution too on the away side for kaons and protons. 

Examining the most central (red symbols, centrality 0 — 10%) and mid central 
(blue symbols, centrality 30 — 40%) events we clearly see a reverse trend up 
to 7 — 10 GeV/c, which is attributed to the region of the minijet production. 
Surprisingly, associated particles on the away side show a stronger effect. 

4 Summary and Outlook 

Monte Carlo studies of two-particle azimuthal correlations have been performed 
with identified trigger particles and identified associated particles in p-p and 
Pb-Pb collisions at mid-rapidity. 

The presented analysis aims to suggest a possible measurement of the charge 
conservation effect via PID-triggered spectra. We think it would be interesting 
to see whether measured data behaves similarly as PYTHIA or HIJING based 
Monte Carlo generated events. The key issues we have found in our analysis: 

• Quantum number conservation can be seen in the near side and we have 
got weak signal on the away side, which suggest the correlations mainly 
arise from parton cascade and hadronization processes in Monte Carlo 
generated data samples. 

• The effect is getting stronger with multiple quantum number conservation 
~ 2 for pions (C), ~ 5 for kaons (C + S), and ~ 10 for proton (C + B). 

• The relative PID-triggered yields defined by Eq. presents increasing 
splitting up to pr < 10 GeV/c and seem saturating at the highest momenta 
at the near side, however the lack of this effect is observed at the away 
side in PYTHIA simulated p-p collisions. 

• In case of HIJING-simulated Pb-Pb collisions this splitting effect is also 
presents, but inverse trend can be seen in the unlike-to-like relative PID- 
triggered yield. This anomalous behavior can be either due to the minijet 
production or high-energy nuclear effect (Cronin, shadowing/anti-shadowing) 
modeled in the HIJING Monte Carlo generator. 

• Both collision energy and multiplicity dependence of the effect were in¬ 
vestigated in p-p in addition to the centrality-dependence simulations 
in central and peripheral Pb-Pb collisions. Here high-momenta or low- 
luminosity events presented larger effect, but this may sensitive to the 
background subtraction (here ZYAM) or the mixture of the soft-hard 
regime in Monte Carlo event generator cocktails. 

In summary our aim is to suggest to perform the presented analysis on real 
data, where either the nature of the correlations or the goodness of the Monte 
Carlo generators can be tested. 
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